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ABSTRACT 

We introduce a new technique that uses galaxy clustering to constrain how satellite galaxies lose 
stellar mass and contribute to the diffuse "intrahalo light" (IHL). We implement two models that 
relate satellite galaxy stellar mass loss to the detailed knowledge of subhalo dark matter mass loss. 
Model 1 assumes that the fractional stellar mass loss of a galaxy, from the time of merging into a 
larger halo until the final redshift, is proportional to the fractional amount of dark matter mass loss 
of the subhalo it lives in. Model 2 accounts for a delay in the time that stellar mass is lost due to the 
fact that the galaxy resides deep in the potential well of the subhalo and the subhalo may experience 
dark matter mass loss for some time before the galaxy is affected. We use these models to predict the 
stellar masses of a population of galaxies and we use abundance matching to predict the clustering of 
several r-band luminosity threshold samples from the Sloan Digital Sky Survey. Abundance matching 
assuming no stellar mass loss (akin to abundance matching at the time of subhalo infall) over-estimates 
the correlation function on small scales (< IMpc), while allowing too much stellar mass loss leads to 
an under-estimate of small-scale clustering. For each luminosity threshold sample, we are thus able 
to constrain the amount of stellar mass loss required to match the observed clustering. We find that 
satellite galaxy stellar mass loss is strongly luminosity dependent, with less luminous satellite galaxies 
experiencing substantially more efficient stellar mass loss than luminous satellites. With constrained 
stellar mass loss models, we can infer the amount of stellar mass that is deposited into the IHL. We 
find that both of our model predictions for the mean amount of IHL as a function of halo mass are 
consistent with current observational measurements. However, our two models predict a different 
amount of scatter in the IHL from halo to halo, with Model 2 being favored by observations. This 
demonstrates that a comparison to IHL measurements provides independent verification of our stellar 
mass loss models, as well as additional constraining power. 

Subject headings: cosmology: theory — dark matter — galaxies: halos — galaxies: structure - 
large-scale structure of universe 



1. INTRODUCTION 

In the concordance ACDM cosmology, galaxies, galaxy 
groups, and galaxy clusters form hierarchically. High- 
density regions condense and virialize, forming bound 
structures known as halos. Halos grow through the 
continual accretion of smaller objects. These accreted 
objects may survive within the virialized region of 
the primary halo as smaller, self-bound, orbiting dark 
matter clumps dubbed "subhalos" or "substructure" 
(iGhigna et all 120001: iKlypin et al.l fl999l: iDiemand et~aH 
2004t iKravtsov et al.ll2004al ). Halos of sufficient mass 



are the natural sites of galaxy formation, with baryons 
cooling and condens i ng towards potential w ell minima 
(jWhite fc Reeslfl978l : iBhimenthal et aLlll986H When a 
halo is accreted by a larger halo, thus becoming a sub- 
halo, the galaxy within it becomes a "satellite" galaxy 
within a group or cluster. Understanding the detailed 
relationship between (satellite) galaxies and (sub) halos 
is a long-standing focus of galaxy formation theory. 

In the hierarchical paradigm, these smaller objects, 
upon merging, become victims of intense tidal fields 
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and interactions within the larger systems in which 
they reside. The dark matter mass associated with 
a subhalo may be rapidly stripped upon infall. This 
stripping acts on the periphery of the subhalo first, 
suggesting that the luminous galaxy, residing in the 
center of the subhalo, may be relatively unharmed. 
After enough time has elapsed, stripping of stars 
may begin to occur as well. These liberated stars 
that are ripped from galax ies are the likely source of 
"intrahalo light " CIHL: e.g.. iGallagher fc Ostrikeil 1 1972 ; 
Merrittl 119831: iBvrd fe Valtonenl ll99Ct iGnedinl l200a 
Murante et"aH 12004); iLin fc Mohrl 12004 iWillman et al l 
2001 IS ommcr- Larsenl 120061: iConrov et all 120071; ? ; 



Purcell et al.l I2007L 120081: ~?T iRudick et al.l I2009L l20ll 

This has been studied in great detail at the scale 
of individual galaxies, o ften categorized as "stellar 
halos" or "tidal streams" (lMorrisqnlll993l: [Sackett et all 



1994t IWetterer fc McGrawl 119961: IMorrison et al.l 



Weil et al.lll997tlChiba feB ccrs 200 01: llvezic et all 



1997) 
20001: 



Leaueux et all 119981: lAbe et all 119991: IMorrison et al" 
20001: lYannv et all 120001: ISiegel et all l2002t llrwin et al 
200a iZibetti fc Ferguson! 120041: iGuhatha kurta et al 
20051: iChapman et all 120061: iKalirai et all 120061: 



McConnachie et all 



Bailin et all l2011h. 



[200l iHood et all [20071: ?: 

galaxy group s (iFeldmeier et all 



20011: ICastro-Rodriguez et all 120031; IWhite et al.l 120031: 
Da Rocha fc Mendes de Oliveira 2005; lAguerri et al.l 
20061 : iFeldmeierl l2006t iDa Rocha et all l2008| ). and 



galaxy clusters (wh ere it is known as the intra-cluster 
light or ICL, e.g.. [Gallagher fc Ostrikeil 119721; iMerrittl 
1983 iMelnick et , al.l 119771: iThuan fc Kormendvl 119771: 
Bvrd fc Valtonenlll990HUson et al.lll991l;lBernstein et all 



19951: ICalcaneo-Rqldan et al.l l200a iGnedinT 12003 



Murante et all 120041: iLin fc Mohrl 12004 IWillman et al 



20041: iMihos et al.H2005t iZibetti et al . 2005; Krick et al 
2006| [S ommcr- Larsenl l200l IConrov et all 120071: 



oy 

Seigar et all 120071: iGonzalez et al.l l2007t"lPierini et all 
20081: iRudick et alll2009Tl2011l: iToledo et alll2011l: ?T~ 



In this paper, we seek to understand the liberation of 
stars from satellite galaxies by connecting stellar mass 
loss 2 to subhalo dark matter mass loss using galaxy clus- 
tering and intrahalo light observations. We employ the 
model for halo substructure introduced in iZentner et all 
(|2005l hereafter Z05) in order to constrain this relation- 
ship. We compare our model predictions to observations 
of the IHL over a large range of host halo mass scales. 
The aim of this paper (Paper I) is to introduce our mod- 
eling framework and its predictive power. In a forthcom- 
ing paper (Paper II) we will extend our analysis to high 
redshift in order to study the assembly of the IHL across 
cosmic time. 

The paper is laid out as follows. In § [5] we discuss 
the motivation for this study. In § [3] we review the Z05 
model for cold dark matter (CDM) substructure. In § 2] 
we describe our models that connect stellar mass loss to 
dark matter mass loss and in §[3] we demonstrate how we 
constrain these models using galaxy clustering. In § [5] 
we show the luminosity dependence of satellite galaxy 
stellar mass loss. In §[7] we use our models to make IHL 
predictions, and compare to observations at low redshift. 
Finally, in § [5] we give a summary of our results and 
discuss directions for future work. 

2. MOTIVATION 

A simple, yet remarkably powerful technique for 
connecting dark matter halo mass to either stellar 
mass or luminosity has emerged in recent years. By 
assuming a monotonic relation between halo mass (or 
maximum circular velocity V ma x = max^GM (< r)/r]) 
and luminosity L (or stellar mass) one can "abun- 
dance match" to make the correspondence between 
dark matter ( sub)halos and an o b servable galaxy 
prope r ty (e.g.. IKravtsov et al.l l2004al: IVale fc Ostrikerl 
2001 iTasItsiomi et all 12004 IVale fc Ostrikerl 120061: 
Conrov et all 120061: IConrov fc Wechslerl l200l 



Mosteretall l2010t iBehroozi et all 120101: iGuo etal l 
20101 : ISimha et al.ll2010l ). For example, this can be done 
by matching the observed number density of galaxies, 
n g , above some luminosity to the number density of 
halos and subhalos, rih, above a certain V ma x, 

n g (> L) = n h (> V max ). (1) 
This yields an implicit relationship between L and V max 
that preserves the observed luminosity function of galax- 
ies. 

2 By "stellar mass loss" , we refer to stars being stripped from a 
galaxy and not gas lost from stars via winds. 



IConrov et all (|2006l ) used this method to assign lu- 
minosities to halos and subhalos in a cosmological N- 
body simulation at several redshifts. They predicted the 
luminosity-dependent clustering of galaxies and found 
that the two-point correlation function (2PCF) of ha- 
los and subhalos matched that of galaxies for a wide 
range of luminosities and redshifts. The authors made 
the physically-motivated choice that Vmax should be the 
maximum circular velocity of subhalos at the time of 
accretio n, VJ^Z, instead of at the time of observation 
(see also iNagai fc Kravtsovll2005l: IVale fc Ostrikerl 120061: 
iBerrier et al.ll2006t ). 

The reasoning behind this choice is as follows. Upon 
merging into a larger host, the Vmax of a subhalo 
will decrease due to mass loss (jHavashi et al.l 120031 : 
IKravtsov et all l"2004al ). Dark matter on the periphery 
of the subhalo will be lost first, because it is less bound 
to the subhalo. On the other hand, the stellar mass of the 
satellite galaxy is concentrated at the center of the sub- 
halo and thus more tightly bound. Tidal stripping can 
significantly alter the surrounding subhalo, but possibly 
leave the galaxy largely unperturbed for some period of 
time (using N-body simulations, ? found that halos need 
to lose ~ 90% of their original mass before tides begin to 
affect the kinematics of stars). Consequently, while Vnax 
of the subhalo decreases, the stellar mass of the galaxy 
may remain unchanged for long periods of time. It fol- 
lows that galaxy observables such as luminosity or stellar 
mass should correlate wit h V^. instead o f the fi nal Vmax 
(|Nagai fc Kravtsovll2005T ) . IConrov et all (|2006l ) lent em- 
pirical support to this picture by showing that the choice 
of V^ax was essential in order to achieve agreement with 
the observed clustering of galaxies. 

It is well known that galaxy groups and clusters are re- 
plete with diffuse stellar material. This material is widely 
thought to be the remains of disrupted satellites (see §[T]). 
However, the choice to associate galaxy luminosity with 
Vnax at the time of satellite accretion is tantamount to 
assuming that no stars become unbound from satellite 
galaxies. In this case, why can using V^™ accurately 
reproduce the observed clustering? 

There is a crucial subtlety to simply using for 
abundance matching. When matching galaxy and sub- 
halo number densities with the L — Vmax relation by 
choosing V^ K , there is an inherent second threshold 
in the final V na x of subhalos. This is a threshold be- 
low which objects are "operationally removed from con- 
sideration" due to the limited resolution of the sim- 
ulation that is used to pe r form the calculation. In 
the study of IConrov et all (j2006f ) . this threshold was 
~ 80km s _1 - the completeness limit of halos in their 
simulation. Therefore, an accreted subhalo will be ar- 
tificially destroyed if it becomes sufficiently small such 
that its structure is not well resolved within the sim- 
ulation. These so-called "orphans" are neglected and 
they unwittingly act to model the stripping of stars. The 
galaxies that would have been associated with these ha- 
los, had they not become unresolved, are effectively re- 
moved from the final galaxy sample, just as if they were 
disrupted. This elimination of subhalos directly affects 
the 2PCF, which is ver y sensitive to subhalo abundance 
(see lWatson et alll2011l ). Including these orphans should 
result in a boost of the small-scale correlation func- 
tion, implying that pure abundance matching using V x 
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max 
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may not accurat ely reproduce the observed clustering on 
small scales (seelKitzbichler fe Whltell2008t iMoster et al.l 
I2010t IWetzel fc Whitell2010D . In fact, we have performed 
a test in which we selected a T^ax > 210km s _1 (corre- 
sponding to ~ L* galaxies and brighter) threshold with 
two, secondary final V max thresholds, V max > 20km s _1 
and Vmax > 80km s — 1 , to mimic this resolution limit 
effect. Increasing this threshold from 20 to 80km s _1 
resulted in a ~ 20% decrease in the 2PCF 3 at scales 
less than IMpc. We speculate that below this resolution 
limit is where much of the stripped stellar material may 
originate. 

The above discussions bring us to the motivation of 
this work. Including stellar mass loss in the abundance 
matching technique is necessary to describe the observed 
clustering of galaxies in a manner that does not include 
implicit selections. Observed galaxy clustering therefore 
has the potential to constrain the typical amount of stel- 
lar mass loss from a subhalo. In this work, we develop 
models that relate satellite galaxy stellar mass loss to 
subhalo dark matter mass loss. We use these models, to- 
gether with the abundance matching technique, to pre- 
dict the observed 2PCF (see § 0] and § [5] for details). 
Comparing to measurements allows us to constrain our 
stellar mass loss models. This requires a detailed under- 
standing of the evolution of subhalos within hosts. We 
use the Z05 analytic model for halo substructure, which 
is not subject to any intrinsic resolution effects. The 
model is capable of tracking surviving subhalos down to 
V max "C 80 km s _1 , so predictions are not affected by the 
"orphan" population of galaxies. Our models yield stellar 
mass loss histories from satellite galaxies, so we can make 
predictions for the amount of intrahalo light at varying 
scales and compare to observations. This provides an 
independent check of our stellar mass loss constraints. 
This modeling can also be done at varying redshifts in 
order to study the evolution of stellar mass loss and the 
assembly of IHL over time. Ultimately, this investigation 
could enable a single model to track the stellar-to-halo 
mass relation as a function of time such that abundance 
matching at multiple redshifts would be unnecessary. 

3. THE SUBHALO EVOLUTION MODEL 

We construct models in which the stellar mass loss 
of a satellite galaxy is linked to the dark matter mass 
loss of its subhalo. Consequently, we require a de- 
tailed understanding of the evolution of subhalos within 
ho sts. For this we us e the Z05 model, which is based 
on IZentner fc Bullockl (120031) and is s i milar t o the in- 
dep endent models of [Taylor fe Babul (|2004 l2005al jbl) 
and iPeharrubia fe Benson! (|2005l ) , and snares many fea- 
tures with o ther approximate t r eatments of subhalo 
populations (|Qguri fc Leel 120041: Ivan den Bosch et all 
120051 iFaltenbacher fc Mathe ws 2005: iPurcell et alll2007f 
IGiocoli et all l2(M 12001 iGan et all \2mR ?). The 
Z05 model produces subhalo mass functions, occupation 
statistics, and radial distributions within hosts that are 
in good agreement with a number of high-resolution N- 
body simulations (Z05, and the recent comparison in 
iKoushiappas et aHl2010T) . 

The analytic model proceeds as follows. For a host 

3 We explain the tools we use to calculate correlation functions 
in §[5] 



halo at a given redshift z, and mass M, we generate a 
halo merger tree using the mass-cons erving implemen- 
tation of the excursion se t fo rmalism dBond et al.|[l991t 
lLacev fc CoIill993lll994l: see IZentnerll2007l for a review) 
developed bv ISomerville fc Kolattl (|1999f >. This pro- 
vides the entire history of all halos that merged to form 
the fin al, host halo. We assign all halos iNavarro et all 
(| 19971 hereafter NFW) density profiles with concentra- 
tio ns determined by th eir merger histories according 
to IWechsler et~aT1 (|2002l ). When a halo merges into a 
larger halo, it becomes a subhalo and is assigned ini- 
tial orbital parameters drawn from distri butions mea- 
sured in N-body simulations (Z05; also see lBensonll2005l 
for similar formalisms). We then integrate each sub- 
halo orbit within the gravitational field of the host, 
subjecting the subhalos to orbital decay via dynami- 
cal friction and mass loss through tides and interac- 
tions. We esti mate dynamical frictio n with an up- 
dated form of the [Chandrasckhar (1943) approxim ation 
([Hashimoto et al.l 120031 : IZentner fcBullockl I2003T ). and 
allow for dark matter mass loss beyond the tidal radius 
on a timescale comparable to the local dynamical time. 
Finally, we account for internal heating so that scaling 
relations de scribing the interna l structures of subhalos 
are obeyed (lHayashi et al.ll2003l : iKazantzidis et aLll2004l : 
iKravtsov et al.ll2004bl ). We refer the reader to Z05 for 
specific details of these ingredients. 

The lack of a resolution limit in the Z05 model allows 
us to track subhalos regardless of how much mass they 
have lost. Therefore, we do not lose subhalos due to 
mass loss. However, we do consider a subhalo to be ef- 
fectively "destroyed" if its orbital apocenter sinks to less 
than r apo < 5 kpc distance from the center of its host 
halo. This choice is physically motivated because the 
galaxy within such a subhalo would likely have merged 
with the central galaxy, or at least be observationally in- 
distinguishable from it. This criterion thus models the 
cannibalism of satellite galaxies by central galaxies. We 
show in § [7] that inclusion of these cannibalized subha- 
los in our predictions for intrahalo light is negligible, so 
this criterion actually has little effect on our modeling re- 
sults. In practice, subhalos lose the vast majority of their 
mass prior to achieving such small apocenters, so this cut 
serves only to terminate the integration of that particu- 
lar orbit. In the end, we amass catalogs of all surviving 
subhalos in the final host halo at the time of "observa- 
tion" and we know exactly how much dark matter has 
been lost from each subhalo. A halo that merges into 
a larger host may contain subhalos of its own. These 
"subs-of-subs" of sufficiently high mass to host an ob- 
servable galaxy are only abundant inside very large host 
masses. They are present in our model, but rare. 

To properly sample the distribution of halos in the uni- 
verse, we compute subhalo populations for a grid of host 
halo masses in the range 11 < log(Af/ft, _1 Af Q ) < 15 (in 
steps of 0.1). To account for statistical variation among 
halos and subhalos, we perform 500 statistical realiza- 
tions of the subhalo population at each host mass. The 
result is 500 host halos along with their subhalos at each 
of 41 distinct masses, giving a total of 20,500 distinct sub- 
halo populations. The model predicts the amount of dark 
matter lost from each subhalo as it orbits in the tidal field 
of its host halo. Figure [1] shows the fraction of subhalo 
mass lost as a function of host halo mass, in bins of accre- 
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Fig. 1. — The average fraction of dark matter lost from subhalos, /dm loss = (-^acc ~~ ^fi™)/^acc> as a function of host halo mass for 
three bins of subhalo accretion epoch ,z a cc (shown by the three point types), and three bins of subhalo mass (at accretion; shown by the 
three panels), according to the Z05 semi-analytic model. Each point represents an average over all subhalos in a given bin of accretion epoch 
and subhalo mass, from 500 model realizations of a specific host mass, and errorbars show the lcr scatter. The points representing different 
accretion epoch bins are slightly staggered for clarity. Accretion time is the dominant factor that determines /dm loss; with subhalos that 
have merged earlier having more time to be stripped of their dark matter. 



tion epoch z acc , and subhalo mass (at accretion). This 
fraction is defined as /dm loss = (M a d ™ - M£™)/M£™, 
where Mf 1 ™ is the subhalo mass at accretion and Mi?" 1 

<acc nil 

is the final subhalo mass. Each panel of the figure 
represents a bin of subhalo mass and the three sets of 
points in each panel represent bins of accretion epoch: 
Zacc = 0.25 — 0.5 (green circles), z acc = 1.0 — 1.25 (red X 
symbols), and z acc = 1.75 — 2.0 (blue diamonds). Each 
point thus represents an average over all subhalos in a 
bin of accretion epoch and subhalo mass, from the 500 
model realizations of a specific host halo mass, and er- 
rorbars show the lcr scatter. It is clear that accretion 
time is crucial towards determining /dm loss, with sub- 
halos that have merged earlier having more time to be 
stripped of their dark matter. Moreover, the scatter in 
/dm loss shrinks for earlier accretion times. At fixed ac- 
cretion epoch, the average value of /dm loss is remarkably 
constant. However, there is a slight increase as we move 
to lower host halo mass at fixed subhalo mass, or as we 
move to higher subhalo mass at fixed host mass. This is 
a result of dynamical friction playing a stronger role for 
subhalos that are considerable in size compared to their 
host halos. The larger a subhalo relative to its host, the 
more rapidly its orbit will decay and it will sink within 
the host potential. The tidal field of the host is stronger 
towards the host center, and this induces more dark mat- 
ter mass loss. 

4. MODELS FOR SATELLITE GALAXY STELLAR MASS 

LOSS 

We now present the models we use to describe satel- 
lite galaxy stellar mass loss. The exact way in which 
stars are stripped from the subhalo they inhabit is likely 
quite complicated. However, we can make some simple, 
physically-motivated approximations aimed at capturing 
the gross, relevant behavior. In this sense, we pursue the 
question of the evolution of stellar mass using a philoso- 
phy similar to that which underlies abundance matching. 
We aim to make a set of minimal, yet effective assump- 



tions that serve to distill the enormous amount of infor- 
mation contained in survey data. Indeed, we aim in part 
to extend the abundance matching techniques by making 
the lower threshold for stellar mass explicit, rather than 
implicit. 

We consider two models in which we relate the amount 
of satellite galaxy stellar mass lost to the correspond- 
ing amount of subhalo dark matter lost. Combined with 
the Z05 model that makes detailed predictions for dark 
matter mass loss, these models can predict the stellar 
mass loss for any given halo. It is important to note 
that the Z05 prescription for mass loss does not include 
the effects that baryons c an have on dark matter . Adi- 
abatic Contraction (AC) (jBlumenthal et aD ll986: ?; ?) 
may increase the central density as the gas condenses 
and sinks to to the cente r of the dark matter potential 
well (jDiemand et al.ll2004t ?; ?; ?). Conversely, processes 
during halo formation, such as gravitational heating from 
merger events, can counteract AC (e.g., ?) and, in fact, 
the central density may decrease through baryonic feed- 
back (?). These effects can impact the survival of satel- 
lite galaxies (see ?), however consideration of these de- 
tailed effects in our modeling is beyond the objective of 
this paper. 

4.1. Model 1 

Our first model sets the fraction of stellar mass that is 
lost from a galaxy to a fixed proportion of the fraction 
of dark matter that is lost from its subhalo. The model 
works as follows. Any halo of sufficiently large dark mat- 
ter mass will have some stellar mass associated with the 
galaxy it contains at the time of accretion, M* cc . This 
stellar mass will be some fraction of the mass in dark 
matter, M~£. After the halo merges into a larger halo, 
becoming a subhalo, it orbits within the host halo po- 
tential and loses mass. At the time of observation, the 
subhalo has a smaller mass, MgJ 1 . We relate the frac- 
tion of stellar mass that is lost during this time to the 
fraction of dark matter mass that is lost through a single 
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Fig. 2. — Schematic of how satellite galaxy stellar mass loss occurs in Model 1 and Model 2. From the time the subhalo merges into the 



host halo, t n 



until the final rcdshift under consideration, ifi na i, our semi-analytic subhalo model predicts the amount of dark matter 



mass loss, AA/ dm = M5* 



-Mi m . 

fin 



Model 1 assumes that the amount of stellar mass lost, AM*, is proportional to AM™. Model 2 defines 



a radius, r', such that stellar mass is only lost if dark matter is lost within that radius. If this is the case, then AM* is only proportional 
to the amount of dark matter lost inside of this radius. See §[4] for detailed descriptions of the models. 
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The left-hand side of Figure [5] is a cartoon schematic of 
how stellar mass loss occurs in Model 1. From the time 
the subhalo merges into the host halo, i mer ge, until the 
final redshift under consideration, t nna i, the Z05 model 
(see §[3] for model details) predicts the amount of subhalo 
dark matter mass loss AM dm with the stellar mass loss 
AM* being related to the dark matter mass loss via e. 
Therefore, we are left with a simple parametric equation 
governed by a single free parameter. For example, if 
e = 0.5, then a subhalo that loses 50% of its dark matter 
will lose 25% of its stellar mass. As we mentioned above, 
stellar mass loss should be less efficient than dark matter 
mass loss, so we should generally expect e < 1. 

In order to compute the final stellar mass of a galaxy 
via Eqn. 2J we need to know the stellar mass of a satel- 
lite galaxy as a function of the subhalo mass at the 



time of accretion. The stellar mass of a galaxy de- 
pends on many physical processes, including mergers, 
gas cooling, star formation, feedback from supernovae, 
feedback from active galactic nuclei, making ah initio 
predictions highly non-trivial. However, there are sev- 
eral empirical methods for obtaining the stellar-to-halo 
mass relation (SHMR), which have recen t ly appeared 
in th e literature Ce.g.. lYang et all l2003t IWang et all 



2006 : iConrqv fc Wechsleri 120091: iWang fc Jind 



Moster et alJl2010HGuo et al.lteOlOHBehroozi et all 
Neistein et all 12011( 1 



We assign stellar mas ses to ha- 
los with the SHMR of iBehroozi et all (|2010l hereafter 
B10). Specifically, we employ the relation used in Eqn. 21 
of B10, with the mean parameter values given by the 
H = k = model in their Table 2. Therefore, for ev- 
ery halo and subhalo in our catalog, we can use the 
halo mass at the time of accretion to assign a stellar 
mass to the galaxy hosted by the halo. The top panel 
of Figure [3] shows how stellar masses are related to our 
host halos and subhalos with the B10 SHMR. The black 
curve shows the mean of the SHMR and grey points il- 
lustrate the assumed 0.15 dex scatter given in B10. We 
show one million points randomly drawn from the full 
distribution of halos and subhalos over the host halo 
mass range 11 < log (M /h~ 1 Mo) < 15. Stellar mass 
rapidly increases as a function of halo mass at low masses 
before turning over and becoming shallower at higher 
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Fig. 3. — The stellar-to-halo mass relatio n (SHMR) for our host 
halos and subhalos with the adopted Bchroozi ct al. (2010) SHMR. 
The black curve shows the mean of the SHMR and grey points il- 
lustrate the assumed 0.15 dex scatter by sampling 1 million halos 
and subhalos from all of the model realizations. The bottom panel 
is the same as the top with the y-axis divided by halo mass to high- 
light the characteristic mass (~ 10 12 h~ 1 Mq) where star formation 
is most efficient. 



host masses. The bottom panel is the same as the top 
with the y-axis divided by halo mass in order to high- 
light the characteristic mass where this turnover occurs, 
-Whaio ~ 10 12 /i _1 A/q. This characteristic mass is the 
halo mass at which star formation is most efficient. 

For most applications, SHMRs are developed by link- 
ing the halo mass function to the galaxy stellar mass 
function (SMF) through abundance matching. This is 
typically done using the subhalo mass at accretion which, 
by definition, presumes that no stellar stripping occurs. 
However, even if stellar stripping occurs (i.e., if e is non- 
zero) we can still use a SHMR relation that assigns stel- 
lar mass at the time the subhalo merges, due to the fact 
that the SMF is strongly dominated by central galax- 
ies. We emphasize that while there are many SHMRs 
in the literature, we are not very sensitive to the partic- 
ular choice of SHMR. As we will discuss in § we as- 
sign stellar masses to halos and subhalos and then rank 
them in stellar mass in order to find the stellar mass cut- 
off (Mg n ) that matches the observed number density of 
a given galaxy luminosity threshold sample. Therefore, 
two SHMRs that yield the same rank order for the halos 
and subhalos will be indistinguishable from each other, 
even though the values of MS will be different (as well 
as the mass-to- light ratios). All published SHMRs are 
monotonically increasing functions, so this behavior is 
general. On the other hand, the scatter in the SHMR is 
important, as this will change the rank order. We use the 
BIO scatter of 0.15 dex throughout our analysis, though 
we test the effect of changing the scatter on our results 
in § 15.11 We also note that we do not consider the fact 
that the stellar mass of a galaxy may actually increase 



for some time after merging and becoming a satellite, as 
it has been recently shown that star forma tion in active 
satell ites may continue for several Gyr (see IWetzel et al.l 
1201 If ). Such detailed modeling is beyond the scope and 
intention of this work. 

In summary, we take the following steps in Model 1 for 
assigning final stellar masses to an ensemble of halos. 

1. We use the Z05 model for subhalo evolution to 
obtain a list of halos and subhalos for a range of 
host halo masses, and we determine the fractional 
amount of dark matter mass lost from each sub- 
halo. 

2. We use the B10 stellar-to-halo mass relation to as- 
sign stellar masses to host halos and subhalos at 
the time of accretion. 

3. For every subhalo, we compute a final stellar mass 
by relating the fractional amount of stellar mass 
lost to the fractional amount of dark matter mass 
lost through the single free parameter e, as de- 
scribed in Eqn. 0] 

4.2. Model 2 

Model 1 operates under the basic assumption that the 
fractional amount of satellite galaxy stellar mass loss is 
proportional to subhalo dark matter loss through the free 
parameter e. This means that even a small amount of 
dark matter stripping will be accompanied by some stel- 
lar stripping. However, as we argued in § [2l a subhalo 
may lose a considerable amount of dark matter before its 
galaxy is significantly disturbed. For example, subhalos 
in high-resolution simulations lose a significant fraction 
of mass at their outskirts rapidly upon merging into a 
larger system, while their interiors remai n unaltered by 
this mass loss (e.g.. iDiemand et"aL| [2007V Accordingly, 
we consider an alternative model that incorporates a de- 
lay between the initial loss of dark matter mass and stel- 
lar mass loss. 

Our second model states that stellar mass will only be 
lost if a sufficient amount of dark matter is lost first. The 
model works as follows. We first define an approximate 
radius of influence that a galaxy has within its subhalo, 
corresponding to the region within which the gravity due 
to the stellar component is comparable to that from dark 
matter. To estimate this radius, we first assign a stellar 
mass to each subhalo at the time of accretion, M* cc , using 
the B10 formalism. Next we assume, for simplicity, that 
this is a point mass and we do not consider the total 
amount of "cold baryons" or assume a galaxy profile. 
To find the radius of influence, r; n fl, of the stellar mass 
associated with the galaxy, we calculate where the dark 
matter mass enclosed within the subhalo is equal to M* cc . 
We thus integrate the NFW mass profile, 



Mace (< Hnfl) = 4 W 3 



In 1 



Hnfl 



iflA 



1 



nfl/r s 



(5) 

where po and r s are the NFW parameters, and solve for 
r infl by setting M a d ™(< r infl ) = M a * cc . 

The radius finfl gives the rough scale within which mass 
must be lost before significant stellar mass is lost; how- 
ever, we allow for flexibility in this final prescription for 
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stellar mass loss. We assume that a galaxy will start 
losing stellar mass once dark matter is stripped from in- 
side a radius that scales linearly with r m n. We define a 
new radius, r' = ip x r; n fl, where ip is a free parameter 
of order unity. If the subhalo loses so much dark mat- 
ter that its final mass is less than Mf™(< r'), we allow 
stellar stripping to occur. Furthermore, we assume that 
the fraction of stellar mass lost is equal to the fraction 
of dark matter lost after the bound mass of the subhalo 
crosses below this threshold. The final stellar mass of the 
galaxy is thus 



M* 



1 



0, 



M*£(< r') 



(6) 

which is analogous to Eqn. [4] 

The right-hand side of Figure [2] illustrates how Model 2 
works. After accretion, if the subhalo shrinks enough due 
to mass loss such that its mass at the final redshift of in- 
terest is less than Mf™(< r'), then stellar mass will be 
lost (AM*) in proportion to the additional amount of 
dark matter lost. However, if mass loss does not reduce 
the subhalo mass to less than Mf™(< r'), no stellar mass 
is lost. Model 2 incorporates a delay between the strip- 
ping of the outer layers of dark matter and the stripping 
of the stellar material residing in the depths of the sub- 
halo potential well. The single free parameter of Model 2, 
ip, allows us to vary the amount of lag. As was the case 
for e = in Model 1, setting ip = in Model 2 leads to 
zero stellar mass loss in all subhalos. 

In summary, we take the following steps in Model 2 for 
assigning final stellar masses to an ensemble of halos. 

1 . We obtain a population of halos and assign stellar 
masses at the time of accretion as described in the 
first two steps of the Model 1 summary. 

2. For each subhalo, we calculate ri n a using Eqn. [SJ 
We then scale this radius using the free parameter 
ip, to get r' = ip x rinfl. 



3. We calculate the dark matter mass at accretion 
that is enclosed by r', M^(< r'). If insufficient 
dark matter loss has occurred, such that the final 
subhalo mass is greater than this, then we assume 
that no stellar mass loss has taken place. Other- 
wise, we estimate a final stellar mass using Eqn. [6l 

It is important to note the effect of allowing the stel- 
lar and dark matter components to have distinct density 
profiles. Using N-body simulations, ?) studied the tidal 
evolution of satellite galaxies embedded in dark matter 
halos. They adopted the observationally-motivated as- 
sumption that stellar components have cored (King) pro- 
files, while dark matter follow cuspy NFW models. They 
helped explain the extremely large mass-to-light (M/L) 
ratios of ultrafaint MW dwarfs by showing that the cuspy 
nature of the NFW profile keeps the central dark mat- 
ter component more tightly bound than the stars. For 
our Model 2, we have taken e = 1 for simplicity in order 
to illustrate the basic features of a simple model as well 
as the manner in which it can be constrained. Taking 
e = 1 in this model allows for stellar mass loss with a 
particular time dependence. Specifically, M/L ratios can 
evolve only when mass loss does not occur interior to the 



threshold radius r' . Any additional mass loss interior to 
r' keeps the M/L ratios fixed. There are two considera- 
tions to this end. First, allowing e to take on a value other 
than 1 would give a time-dependent M/L ratio even for 
mass loss interior to r' . However, our goal is to present a 
simple model and illustrate constraints on such a model. 
A second parameter would lead to inevitable degenera- 
cies among parameters and introduce complexities that 
would detract from the broader point of the current pa- 
per, and we show below that a single-parameter model 
adequately describes various observational data. Second, 
in even greater detail, we could consider a time-varying 
e = e(t), to describe a very broad range of variations in 
M/L ratios with time. Though not currently necessary to 
describe the data we consider, such a model may be more 
physically appropriate, and the time variation of e could 
be derived, at least approximately, from specific assumed 
stellar mass profiles. In follow-up work, we are building 
on the complexity of the models we present here, though 
the details of these broader models are well beyond the 
scope and intent of this paper. 

5. CONSTRAINING SATELLITE GALAXY STELLAR MASS 
LOSS USING GALAXY CLUSTERING 

We now turn to constraints on the relationship be- 
tween stellar mass and halo mass imposed by the cluster- 
ing of galaxies. Our two stellar mass loss models spec- 
ify the final stellar mass in any halo or subhalo, given 
a value for the parameter e (for Model 1) or ip (for 
Model 2). With stellar masses assigned to all halos and 
subhalos, we can predict a ha l o occ u pation distribution 
(HO D: e.g.. iPeacock fc Smith! 120001: IScoccimarro et al.l 



ifM iBerlind fc Weinberg! 120021: iCoorav fc Shethl 12002ft 
for any given stellar mass threshold. In particular, we 
compute the mean number of galaxies as a function of 
host halo mass, (N)m- We calculate the mean over the 
500 host halo realizations at each host mass. We use 
this function to compute the number density of galax- 
ies by weighting the host halo abundances by (N)m and 
integrating over all halo masses, 



f°° dM-^-(N) 
Jo dM x 1 



M ■ 



(7) 



We adopt the IWarren et all (|2006T) halo mass function, 
dn/dM, in this calculation, though our results are not 
sensitive to the specific choice of mass function. In this 
way, we find the stellar mass threshold that yields a 
galaxy number density equal to that of the observed sam- 
ple with which we aim to compare. 

We compare our model pre dictions to four of the Sloan 
Digital Sky Survey (SDSS: lYork et al.l f2000h luminos- 
ity th reshold galaxy samples measured by IZehavi et al.1 
(|2011[ ). Specifically, we consider volume-limited sam- 
ples with r-band absolute magnitude thresholds of M r < 
— 18, —19, —20, and —21 (nu mber densities for these 
samples are listed in Table 2 of lZehavi et"alll20111 ). We 
do not consider the fact that the observed relationship 
between stellar mass and luminosity is not one-to-one 
(i.e., galaxies of a fixed luminosity can have different stel- 
lar masses which, for instan ce, is manifested as sc atter in 
the Tully-Fisher relation - iBell fc de Jong] 1200 1[ ). How- 
ever, including this scatter should be similar to increasing 
the scatter in the B10 relation between dark matter mass 
and M* cc , which has a very small effect, as we show in 



Once we have determined the appropriate stellar mass 
threshold, we use the resulting HOD to calculate the 
predicted clustering of the model. We place galaxies 
within the host halos in a N-body simulation to mea- 
sure the 2PCF. We use a single realization of the "Con- 
suelo" simulation (with a box size of 420/i~ 1 Mpc and 
1400 3 particles), which is part of the LasDamas suite 
of simulations (McBride et al., in prep.). The cosmol- 
ogy assumed in our semi-analytic subhalo model is set to 
match that of the LasDamas simulations 4 and is si milar 
to the recent WMAP7 values (jKomatsu et alj|2011f) . We 
use t he spherical over-den sity (SO) halo finder "Rock- 
star" (jBehroozi et al.ll2~011h to identify halos. We use an 
SO halo finder because the SO algorithm mimics the as- 
sumptions made in our semi-analytic model. Both the 
Z05 model and the halo finder define v irial masses based 
on th e virial threshold definition of iBrvan fc Normanl 
(1998). We note, however, that our results are not sen- 
sitive to the choice of halo finder. We have repeated our 
analysis using a Friends-of-Friends halo finder and de- 
rived similar results. We populate host halos with galax- 
ies according to our (N)m model predictions. 

In order to compare to data, we convert our real- 
space correlation functions, £(r), to projected correla- 
tion functions, w„ (r D ), by integrating alon g lines of sight 
(|Davis fc Peebledll983t IZehavi et alJl2004j ): 

Wp (r p ) = £(y/t*+v*)dy. (8) 

For each lum inosity threshold, w e integrate out to the 
same y max as IZehavi et al.l ([20111 listed as 7r max in their 
Table 2). We predict w p (r p ) for the four luminosity sam- 
ples, given any value of e from Model 1 or ip from Model 2. 
In order to constrain our stellar mass loss models, we 
run a grid of e and ip values in steps of Ae = 0.1 and 
Aip = 0.05. For each parameter value and each luminos- 
ity threshold, we compute y 2 by com paring our model 
prediction to the IZehavi et al.l (|2011l) measurement of 
Wp^p), using only diagonal errors (not the full covari- 
ance matrices). 

Figure 0] shows clustering r esults for Mo d el 1 com- 
pared to measurements from Zeh avi et al.l (|2011| ) for 
our four SDSS luminosity threshold samples: M r < 
— 18, —19, —20 and — 21. To reiterate, Model 1 assumes 
that the fractional amount of satellite galaxy stellar mass 
loss is proportional to the fractional amount of subhalo 
dark matter mass loss with proportionality constant e. 
Adopting e = specifies no stellar mass loss and is equiv- 
alent to using Vmax at the epoch of accretion for abun- 
dance matching. On the other hand, adopting e = 1 
means that stellar mass is lost at the same rate as dark 
matter. For each luminosity sample, we show the e = 
and 1 cases as red solid and dashed lines, respectively. 

At each luminosity, e = predicts small-scale clus- 
tering that is stronger than the data indicate and this 
gets to the essence of using clustering to constrain stellar 
mass loss. This over-prediction can be attributed to the 
fact that no stellar mass loss has occurred for satellite 
galaxies, resulting in too many satellites with high stel- 
lar masses. This leads to enhanced clustering on small 

4 Throughout the paper, we work within the standard, vacuum- 
dominated, cold dark matter (ACDM) cosmological model with 
Q m = 0.25, n A = 0.75, tt h = 0.04, h = 0.7, erg = 0.8, and 
n s = 1.0. 
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Fig. 4. — Clustering predictions for Model 1. Model 1 sets the 
fractional amount of satellite galaxy stellar mass loss to be propor- 
tional to the fractional amount of sub halo dark matter mass loss 
through the free parameter e (see H4,ll for model details and Fig. [2] 
for an illustration). Each panel shows results for a different galaxy 
luminosity threshold. Solid red curves show model predictions for 
e = 0, equivalent to no stellar mass loss. Dashed red curves show 
the e = 1 case, where the stellar mass loss occurs at the same 
rate as subhalo dark matter mass loss. Grey curves show the e 
values that m atch the observed SDSS clustering measurements of 
Zehavi et al. (2011, black points). 

scales (e.g. IWatson et all 1201 lh . To be consistent with 
the data, some stellar mass loss needs to occur. How- 
ever, the e = 1 result demonstrates that too much stellar 
mass loss leads to weaker clustering than the data re- 
quire. The grey curves in Figure |4] correspond to inter- 
mediate amounts of stellar mass loss and represent the e 
with the lowest \ 2 values. 5 Despite the fact that Model 1 
paints a simplified picture of how stars are stripped from 
the galaxies they reside in, we show that it is very effec- 
tive at matching the observed clustering. Moreover, it is 
striking that e decreases with increasing luminosity. We 
will discuss this luminosity dependence in detail in § |6] 

Figure [S] shows similar results for Model 2. This model 
was designed to allow dark matter on the periphery of 
a subhalo to be lost due to the strong tidal field of the 
host, without significantly altering the luminous galaxy 
residing deep in the core of the subhalo. Model 2 mimics 
this "lag" by defining a critical radius within each sub- 
halo and only allowing for stellar mass loss if the subhalo 
loses mass from within that radius. This radius can be 
varied through the parameter ip. The limit ip = corre- 
sponds to no stellar mass loss and is equivalent to abun- 
dance matching using Kn ax at accretion. In Model 1, 
e = 1 is analogous to using the subhalo mass at the fi- 
nal redshift output for abundance matching. In Model 2, 

5 We note that none of our model predictions match the large- 
scale clustering for the M r < —18 sample. This may be at- 
tributed to the finit e volume of the sample, as discussed in § 3.2 of 
IZehavi et all iWTH) . 
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Fig. 5. — Clustering results for Model 2. Model 2 defines a 
radius of influence, r; n g, for a galaxy within a subhalo. If a subhalo 
loses sufficient dark matter mass such that it starts to lose mass 
inside of r; nH , then stellar mass is lost at a rate proportional to 
the subsequent subhalo mass loss. The free parameter ip allows 
rinj to grow or shrink, thus a larger value of ip means that more 
stellar mass loss will occur on average for satellite galaxies. Solid 
red curves show model predictions for ip = 0, equivalent to no 
stellar mass loss, thus the same as e = for Model 1. Dashed red 
curves show the predictions for ip = 2, an arbitrary value chosen to 
represent how too much stellar mass loss will under-predict ui p (r p ). 
Grey curves show the ip values that match the observed clustering 
for the same luminosity threshold samples as in Fig. [4] 



this happens when the critical radius is larger than the 
size of the subhalo at accretion. This occurs at a dif- 
ferent value of ip for each subhalo. Thus, we choose an 
arbitrary value ip = 2 to again illustrate how too much 
stellar mass loss will under-predict clustering. As was 
the case for Model 1, we are able to match the observed 
clustering and we see a strong evolution in ip as a func- 
tion of luminosity. This emphasizes what was found for 
Model 1, that low-luminosity galaxies experience more 
efficient stellar stripping throughout their evolution in a 
host halo than luminous galaxies. 

5.1. The Effect of Scatter in the Stellar-to-Halo Mass 

Relation 

Observations indicate that halos of a given mass can 
host galaxies w ith a range of luminosit i es and stel- 
lar masses (e.g.. [van den Bosch et all 120071 : IZheng et al.l 
120071: lYang et alj|2008t iMore et all 120091) . For this rea- 
son, it is overly restrictive to assume a one-to-one rela- 
tion between halo mass and stellar mass, whereas a rela- 
tionship with some intrinsic scatter is more appropriate. 
This is especially true for abundance matching because 
the presence of scatter changes the rank order of galax- 
ies. We rank our halos and subhalos by stellar mass, so 
the rank order will be affected by the assumed amount of 
scatter. Scatter is typically accounted for by assuming a 
distribution of stellar mass at fixed halo mass, with the 
mean value of stellar mass given by a particular SHMR. 
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Fig. 6. — The effect of the scatter in the stellar mass - halo mass 
relation, on the clustering of a stellar mass threshold galaxy sample. 
Top panel: The grey curve shows the predicted correlation function 
from Model 1 with e = 0.3 and ou r assumed scatter of a = 0.15 
dex, given by Bchroozi et al. (2010). Solid and dashed black curves 
show results for changing the scatter to a = 0.1 dex and 0.2 dex, 
respectively. Each curve is averaged over three realizations of that 
amount of scatter. Bottom panel: The percent deviations in the 
correlation function from the fiducial scatt er of a = 0.15 . At the 
scale of the innermost data point of Zeh avi e t al. ( 20111') that we 
compare our model results to (~ 0.17h _1 Mpc), lowering the scatter 
to c — 0.1 increases w p (r p ) by ~ 8% and increasing the scatter to 
o~ = 0.2 lowers ui p (r p ) by ~ 5%. On larger scales the effect is much 
smaller. 



For example, the B10 scatter that we adopt in this paper 
is a log- normal distribution with a dispersion a = 0.15. 
We implement this scatter by randomly drawing from 
this distribution when we assign stellar masses to halos 
and subhalos at the time of accretion. The results shown 
in Figures |4] and [5] include this scatter. 

We repeat our analysis with different dispersions in 
order to test the sensitivity of our results to uncertainty 
of the scatter in the SHMR. Figure [5] Shows the effect of 
changing the SHMR dispersion for one of the models that 
matches the clustering of M r < —21 galaxies: Model 1 
with e = 0.3. This sample is the most sensitive to scatter 
because the slope of the SHMR is shallowest for luminous 
galaxies and thus a given amount of scatter in stellar 
mass at fixed halo mass translates into a large scatter in 
halo mass at a fixed stellar mass threshold (see Fig. 
The top panel of Figure [6] shows the correlation function 
for three SHMR dispersions: a = 0.1, 0.15, and 0.2. The 
bottom panel depicts percent deviations from the fiducial 
scatter of a = 0.15. This range is conservative, as it is 
larger than the errors in scatter quoted in B10 (especially 
at the low end). The figure demonstrates that changing 
the scatter has a small effect on w p (r p ), in the sense 
that more scatter generally leads to a lower amplitude of 
clusterin g on small sca l es. O ver the range of data points 
from the lZehavi et al.l (|2011[ ) data that we model, there 
is a maximum ~ 8% and ~ 5% shift in w p (r p ) for the 
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Fig. 7. — The luminosity dependence of satellite galaxy stellar mass loss. Left panel: The best-fit Model 1 parameter e as a function of 
absolute r-band magnitude. The decreasing trend of e with increasing luminosity means that low-luminosity satellite galaxies experience 
greater stellar mass loss relative to subhalo dark matter mass loss than luminous galaxies. Right panel: The best-fit Model 2 parameter 
ifi as a function of absolute r-band magnitude. The decreasing trend of ip with increasing luminosity means that low-luminosity satellite 
galaxies have a greater radius of influence (causing more stellar mass loss) than luminous galaxies. In other words, the same qualitative 
luminosity trend is seen in both models. 



a = 0.1 and a = 0.2 cases, respectively. This is not a 
large effect and does not result in a significant change to 
the best-fit values of e or ip. Uncertainties in the scatter 
of the SHMR do not significantly alter our results or our 
primary conclusions. 

6. LUMINOSITY DEPENDENCE OF SATELLITE GALAXY 
STELLAR MASS LOSS 

Figures [4] & [5] indicate that the best-fit values of both e 
and ip vary with luminosity. We highlight these trends in 
Figure [7J The left-hand panel summarizes results from 
Model 1, and demonstrates that the best-fit value of e 
drops from e = 0.8 for M r < —18 galaxies to e = 0.3 
for M r < —21 galaxies. There is a strong luminosity 
dependence associated with satellite galaxy stellar mass 
loss, such that relatively low-luminosity satellite galaxies 
lose more stellar material per unit of dark matter than 
more luminous satellite galaxies. The right-hand panel of 
Figure [7] summarizes results from Model 2, which show a 
consistent trend. The best-fit value of ip drops from ip = 
0.95 for M r < -18 galaxies to ip = 0.35 for M r < -21 
galaxies. Lower values of ip indicate less stellar mass loss 
for a given amount of dark matter mass loss. The two 
models we explore absorb a significant number of subtle 
effects into simple assumptions. Model 1 and Model 2 
differ substantially in detail, yet both analyses indicate 
that the luminosity dependence of galaxy clustering can 
be explained if low-luminosity satellite galaxies lose stars 
more efficiently than luminous galaxies. 

What does this result mean? Presumably it means 
that the stars within luminous galaxies are more tightly 
bound - relative to their surrounding dark matter - than 
the stars in less luminous galaxies. This could be due to 
luminous galaxies having more compact stellar density 
profiles. In fact, it may be that this result can be ex- 
plained by the well-known correlation between galaxy lu- 
minosity and morphology, whereby luminous galaxies are 
more likely t o be ellipticals t han lo w-luminosity galaxies. 
For example, iBlanton et al.l (|2003|) showed that the mean 
Sercic index of galaxies roughly doubles from M r = — 18 



to M r = —21 (looking at their Fig. 9 and converting i- 
band to r-band magnitudes). We must, of course, keep 
in mind that these morphology-luminosity correlations 
apply to all galaxies, of which satellites are only a small 
portion. It may be that these correlations would vanish 
in samples containing only satellite galaxies. 

An ancillary consideration is the relative size of galax- 
ies compared to the typical (sub)halos they occupy. 
For instance, usin g a sample of 140,000 SDSS galaxies, 
IShen et all ()2003l ) showed the size distribution (half-light 
radii) of early- and late-type galaxies as a function of 
r-band luminosity and stellar mass. Examining their 
Figs. 2 & 3, the half-light radius R increases roughly by a 
factor of - 3 from M r < -19.25 to -22.25 for late-type 
galaxies, and by a factor of ~ 5 for early types. However, 
if we consider the typical halo masses that host these 
galaxies (s ee M m - ln values for M r < — 19 and —22 from 
Table 3 of lZehavi et al.ll2011[ h they correspond to virial 
radii of ~ 125/i~ 1 kpc and ~ 750/i _1 kpc, respectively, 
a factor of 6. In other words, at a very approximate 
level, larger (more luminous) galaxies are slightly deeper 
within their host halo potential wells than smaller (less 
luminous) galaxies, and are relatively less susceptible to 
stellar stripping. 

The luminosity trend that we have found demonstrates 
the power of our modeling approach. We use one observa- 
tional measurement, galaxy clustering, to shed light on 
an important physical process that is essentially unob- 
servable: stellar mass loss from satellite galaxies relative 
to dark matter mass loss. In the next section we show 
how our models can be used to predict a different ob- 
servable, the IHL. The IHL is a remnant of stellar mass 
loss, so it may be used to probe the liberated stellar de- 
bris, and thus provide an important cross-check on our 
results. 

7. PREDICTIONS FOR INTRAHALO LIGHT AT VARYING 

SCALES 

If stars are freed from their host galaxies they will pre- 
sumably become part of the IHL. Some of the IHL is 
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Fig. 8. — /ihl+bhg predictions for the two extreme cases of Model 1: e = and 1. /ihl+bhg is the combined amount of stellar mass 
from the brightest halo galaxy (BHG) and the intrahalo light (IHL) divided by the total stellar mass of a system (the BHG, IHL, and the 
stellar mass still associated with surviving subhalos). Thick solid black curves show the mean /ihl+bhg as a function of host halo mass 
and are computed from 500 model realizations of each host mass, while grey curves represent the median. Light and dark purple shaded 
regions represent the 68% and 95% range of th e model distributions. BHG and IHL contributions are shown individually as dashed and 
dotted curves, respectively. Red points show the Gonzalez ct al. (2007) data measurements for comparison. Left panel: e = means that no 
stellar mass loss occurs, resulting in no IHL. As a result, /ihl+bhg i s strongly under-predicted at cluster scales (M^aio ~ 14 — 15h~ 1 Mq) 
relative to the data. Right panel: e = 1 means that stellar mass loss occurs at the same rate as subhalo dark matter mass loss. This 
re sults in an over-predic tion of the IHL. Note: data points are shifted by ~ 40% in order to convert the Afsoo (A cr ; t = 500) masses given 
in lGonz alcz ct al. (200J) to our virial masses (using the definition A v ; r = 377). 



a diffuse background, while some may be composed of 
coherent streams, a prominent example of which is the 
Sagittarius Stream in our own galaxy (e.g.. Ilbata et ail 
l200lHDohm-Palmer et al.ll200H iNewberg et al.ll2002t ). if 
the stars were stripped from a subhalo in the past few 
dynamical times. Observations o f the I HL are notori- 
ously difficult, and lGonzalez et all (|2007t ) argue that the 
proper quantity to measure is the sum of the light from 
the host central galaxy (the "brightest halo galaxy" or 
BHG) and the diffuse intrahalo light component, because 
it is non-trivial to disentangle the two. Our models as- 
sign stellar masses to all halos (both hosts and subhalos) , 
so we can compute the contribution of the BHG to the 
total stellar mass of the system. This combination is also 
less subject to theoretical uncertainty because collisions 
between the BHG and an infalling satellite can disperse 
stellar mass into the IHL. While parsing of the stellar 
mass of the BHG and IHL individually is sensitive to the 
fraction of stellar mass ejected in such collisions, the sum 
is robust to uncertainties in this process. 

We start by investigating the predictions for the e = 
and e = 1 cases from Model 1. These extremes were 
shown in § [5] to over- and under-predict the observed 
clustering, respectively, for each luminosity threshold. In 
Figure [H we show the fraction of the total stellar mass 
of a system (including stellar mass in satellite galaxies), 
that is contained within both the IHL and the BHG, 
/ihl+bhg, as a function of host halo mass, -Mhaio- The 
thick solid black curves in Figure |8] show the mean value 
of /ihl+bhg obtained from our 500 realizations of host 



halos, and grey curves represent the median. The in- 
ner and outer shaded regions represent the 68% and 95% 
ranges of the model distributions. The mean individ- 
ual BHG and IHL contributions are shown distinctly as 
dashed a nd dotted curves, resp ectively. The /ihl+bhg 
data from Gonzale z et al.l (|2007l ) are shown as red circles 
with errorbars, and are the same in both panels. 

For e = there is no stellar mass loss, and thus no con- 
tribution to the total stellar mass of the system from the 
IHL. The stellar mass of the BHG alone provides a decent 
description of the data for halo masses below group mass 
scales, Mhaio J$ 10 13 h~ 1 M Q . While less massive systems 
have been observed to have some detectable IHL (see §[!]), 
it is mainly large group- and cluster-size objects that are 
known to have significant IHL components. It is there- 
fore not surprising that the e = model underestimates 
/ihl+bhg on these scales. 

At the other extreme, the e = 1 model assumes that 
stellar mass is lost with the same efficiency as dark mat- 
ter. We argued in § 2] that this over-estimates stellar 
mass loss, so it is not surprising that choosing e = 1 
leads to an overestimate of /ihl+bhg for systems greater 
than large groups (M ha io > 10 135 h~ 1 M Q ). At lower 
masses, typical infalling satellites have halo masses be- 
low M^™ < 10 11 H~ 1 Mq and bring with them only a 
small amount of stellar mass (B10). This is the regime of 
the SHMR where stellar mass rapidly decreases with de- 
creasing halo mass (see Fig. [3]). Data are not yet precise 
enough to be sensitive to the stellar mass carried into the 
systems by such small, infalling subhalos. Consequently, 
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Fig. 9. — /ihl+bhg predictions for our best-fit Model 1 (left panel) and Model 2 (right panel). We compute each subhalo's contribution 
to the IHL according to its initial stellar mass by using the appropriate best-fit values of e or tp. For example, a subhalo with a stellar mass 
that corresponds to a galaxy of absolute r-band magnitude equal to -18.5 loses stellar mass according to e = 0.8 (in Model 1), whereas a 
subhalo with a stellar mass corresponding to Mr = —20.5 loses stellar mass according to e = 0.5. Shaded regions, curves, and points are as 
in Fig. [8] Our stellar mass loss models predict that the mean IHL fraction in halos increases from a few percent at the scale of individual 
galaxies (M^aio ~ 10 12 /i~ 1 M ( 7)) to ~ 20 — 25% at cluster scales (Af^alo ~ 10 14 — 10 15 h~ 'Mq). Both models succeed in matching the mean 
trend observed by Gonzalez ct al. (2007). However, Model 1 seems to under-predict the scatter in /ihl+bhg > whereas Model 2 predicts a 
scatter that is more consistent with the data measurements. We thus find that this more physically- motivated model for satellite galaxy 
stellar mass loss is favored over Model 1. 



/ihl+bhg ~ /bhg and can be adequately described in 
low-mass systems at both the e = 1 and e = extremes. 

We note that in our model predictions for the IHL, 
we only consider the stellar mass stripped from satellite 
galaxies that have stellar masses equivalent to M r < — 18 
galaxies and brighter. We ignore stellar mass loss from 
less luminous (less massive) satellites in the /ihl+bhg 
calculation. However, the stellar mass locked up in these 
lower mass "uncounted" subhalos only amounts to ~ 3% 
of the total stellar mass within a given host halo, so ne- 
glecting them has a negligible effect on our /ihl+bhg 
model predictions. We also note that in Figure [8] we 
have shifted the data points to higher masses by ~ 40% 
in order to conver t the M 500 (A cr j t = 500) masses given 
in lGonzalez et all ()2007l ) to our virial masses (using the 
definition A v j r = 377). We assumed an NFW profile to 
perform this conversion. 

Having shown that these two extreme cases are not 
able to describe the data accurately, we turn to the 
predictions of our models. To make /ihl+bhg predic- 
tions at a given host halo mass scale we need to assign 
each (sub)halo an e value for Model 1 and a ip value for 
Model 2. Let us examine just the Model 1 case for sim- 
plicity. When we found the best-fit value of e for each 
luminosity threshold sample in § [5l we also determined a 
final stellar mass threshold for the sample through abun- 
dance matching. For example, the M r < —20 sample 
was found to correspond to halos and subhalos with fi- 
nal stellar masses greater than some value M^ n 20 , where 
final stellar masses are predicted using e = 0.5. The 
M r < —21 sample corresponds to stellar masses greater 



than Mg n 21 , where final stellar masses are predicted us- 
ing e = 0.3, and so on. We assign each subhalo a value of 
e using these stellar mass thresholds. For example, if a 
subhalo has a stellar mass at accretion that falls between 
Mg n 20 and M| n 21 , we assign it e = 0.5. If, instead, it has 
a mass between M| n lg and M| n 20 , we assign it e = 0.7. 
It is awkward to use the stellar mass at accretion to as- 
sign subhalos to ranges in Mg n , but we cannot estimate 
final stellar masses without having a value of e in the first 
place, so we are stuck with this approximate method. 
We follow the same process for Model 2 by choosing the 
appropriate ip values. We emphasize that these e and 
ijj values are not "tuned" to agree with the /ihl+bhg 
observations, rather are genuine predictions that result 
from parameter values found from matching to clustering 
data. 

Figure M shows the Model 1 (left panel) and Model 2 
(right panel) /ihl+bhg predictions, as well as the in- 
dividual IHL and BHG contributions. While /ihl+bhg 
rapidly decreases with increasing host halo mass, the con- 
tribution to /ihl+bhg from the IHL increases. This sup- 
ports a picture in which galaxy-mass halos have very 
few luminous satellite galaxies and little IHL, so the 
stellar masses of these systems are dominated by the 
BHGs. Galaxy size halos accrete the majority of their 
mass from s mall subhalos that h ave very high mass-to- 
light ratios (jPurcell et al.l 12007). Therefore, shredded 
satellites deposit very little stellar mass into the IHL. 
The number of large satellites increases with host halo 
mass. The larger and more common subhalos within 
larger host halos may form stars more efficiently (have 
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Fig. 10. — Histograms of the number of satellites that contribute 
75% of the total IHL to a host halo of mass 10 14 7i _1 Ar e for 500 
realizations of Model 1 (solid histogram) and Model 2 (dotted his- 
togram). It is clear that in Model 2 the same amount of IHL 
comes from a much smaller number of satellites that have each 
been stripped of more stars than in Model 1. Since fewer satellites 
contribute to the IHL in Model 2, the Poisson scatter in this model 
is larger and is evident in the /ihl+bhg scatter in Fig. |91 



lower mass-to-light ratios) and provide a source for devel- 
oping significant IHL. Both models find that host masses 
of Mhost/fc^MQ = 10 12 ,10 13 ,10 14 , and 10 15 , contain 
mean IHL fractions of f mh ~ 2%, 10%, 20%, and 25%, 
respectively. This is in good agreement with several pre- 
vious studies, as we discuss in §[5] The IHL thus provides 
an independent check of our stellar mass loss results. 

Figure [S] contains a comparison betwee n our model 
predictions and the iGonzalez et al.l ()2007f ) data. Both 
models predict mean /ihl+bhg fractions that are consis- 
tent with the observations. 6 Interestingly, we find that 
the scatter predicted by Model 1 is substantially smaller 
than it is for Model 2 and it appears to be i n consis - 
tent with the scatter found by Gonzale z et al.l (|2007l) . 
On the other hand, Model 2 is more successful in repro- 
ducing the observed scatter in /ihl+bhg- To quantify 
this, we calcula t e the probability that the scatter in the 
IGonzalez et al.l (|2007l ) data points could result from the 
model distributions. We do this as follows. First, we re- 
strict the comparison to the high mass regime (the last 
8 data points) because we wish to investigate the scat- 
ter at fixed mass (if we include lower mass points, the 
measured scatter will be affected by the mean trend of 

6 W e note that this discussion assumes that the Gonzalez et al. 
(2007) data are accurate, whereas they may be systematically 
overestimated due to an assumed constant mass-to-light ratio 
IILeauthaud et aLllMTH ). 



/ihl+bhg vs. mass). We next create a model data set 
by choosing a singl e model realization of a host halo (out 
of the 500) for each IGonzalez et a l. (2007) data point, at 
the appropriate halo masses. This leads to 8 model val- 
ues of /ihl+bhg- We then add a Gaussi an random error 
to the se values, using the corresponding Gonzale z et al.l 
(2007) errorbars. To fully sample the distribution of data 
sets that could arise from our models, we repeat this 
process 10,000 times. Finally, we measure the weighted 
standard deviation of /ihl+bhg values for the observed 
data, as well as each of the 10,000 model data sets. We 
find that in Model 1, only 0.8% of the model data sets 
have a scatter that exceeds that of the Gonzalez et al.l 
(|2007l) data, whereas in Model 2, 25% of the model data 
sets have a larger scatter. 

This statistical test indicates that Model 1 should be 
rejected (at the 2.7a level) in favor of Model 2, an indi- 
cation that is not entirely surprising given that Model 2 
accounts for delayed stellar mass loss that is to be ex- 
pected on physical grounds. Model 1 assumes that stel- 
lar mass loss is simply proportional to subhalo mass loss, 
so the scatter predicted by this model must be primarily 
caused by the scatter in the total amount of dark matter 
mass loss, with a secondary amount of scatter coming 
from the SHMR. However, the "all or nothing" action 
provided by the radius of influence in Model 2 some- 
how amplifies this scatter. Upon closer investigation, we 
find that the difference in scatter between the two models 
arises from simple Poisson statistics. Figure [TU] shows the 
distribution of the number of satellites that contribute 
75% of the IHL to a host halo of mass 10 14 h~ 1 MQ, mea- 
sured from 500 realizations for both Model 1 (solid his- 
togram) and Model 2 (dotted histogram) . It is clear that 
in Model 2 the same amount of IHL comes from a much 
smaller number of satellites that have each been stripped 
of more stars than in Model 1. This is a direct result of 
the nature of Model 2: most subhalos do not lose stars, 
but if they do, it is at a high rate. Since fewer satellites 
contribute to the IHL in Model 2, the Poisson scatter 
in this model is larger and is manifested in the greater 
/ihl+bhg scatter seen in Fig. M 

8. SUMMARY & DISCUSSION 

While galaxy formation is complex, we have shown in 
this paper that an understanding of the fates of satel- 
lite galaxies provides key insight into the galaxy for- 
mation puzzle. Satellite galaxies live extremely tumul- 
tuous lives. Their spatial clustering can shed light on 
how they lose stellar mass and contribute to building 
the IHL over a large range of host halo mass scales. 
Part of our motivation for using galaxy clustering to 
constrain stellar mass loss stems from possible short- 
comings in the com monly-employed abund ance match- 
ing technique (e.g., IWetzel fc White! 120101: ?, Reddick 
et al., in prep.), which suggest natural generalizations 
with the gross form of the particular models we study. 
The standard method of abundance matching involves 
mapping galaxies to halos by assuming a relationship 
between halo (and subhalo) mass (or circular velocity) 
and a galaxy property (luminosity or stellar mass) at a 
particular time. A galaxy luminosity assignment accord- 
ing to halo mass at the time the halo was first accreted 
onto a larger system provides a very usef ul description 
of known galaxy clustering properties fe.g.. lConrov et all 
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2006). However, this type of assignment implies that no 
stellar mass loss will occur. This is contrary to many 
aspects of galaxy, group, and cluster evolution, including 
the prevailing paradigm in which the IHL is produced 
by stellar mass liberated from infalling g alaxies by in- 
teractions in the group environment (e.g.. IConrov et al.1 
l2007tlPurcell et alj2007t ?). Abundance matching is also 
sensitive to the resolution of the cosmological simulation 
on which it is implemented. Subhalos are artificially de- 
stroyed by falling below the numerical resolution limit 
and immediately removed from the galaxy population ir- 
respective of their circular velocities at their accretion 
times. This can lead to differing predictions for the 
small-scale clustering of galaxies of a particular luminos- 
ity among simulations with a fixed mapping of halo size 
onto luminosity. 

These considerations persuaded us to explore simple 
models that associate stellar masses to subhalo masses 
and constrain satellite galaxy stellar mass loss by match- 
ing to clustering data. We implemented two distinct 
models. Model 1 related the fractional amount of stel- 
lar mass loss to that of the subhalo dark matter mass 
loss from the time of accretion until the final redshift 
through a single free parameter e. Model 2 was intro- 
duced to mimic a "lag" in stellar mass loss wherein the 
satellite galaxy at the core of the subhalo experiences no 
stellar mass loss until the subhalo has been stripped of 
sufficient dark matter on its periphery. Stellar mass loss 
in this model only occurs once dark matter has been lost 
interior to a radius of influence of the satellite galaxy. 
We introduced the free parameter ip that dictated how 
the radius of influence could expand or shrink allowing 
for more or less stellar mass loss. We calculated the val- 
ues of e and ip that matched the observed clustering of 
SDSS galaxies over a large range of luminosity threshold 
samples. Matching to the observed clustering at each 
luminosity threshold directly informs us of the amount 
of stellar mass lost from satellite galaxies and, thus, will 
comprise the IHL. 

This procedure enables predictions for the amount 
of IHL or IHL+BHG compared to the total stellar 
mass of a host system of a given mass (/ihl and 
/ihl+bhg, respectively). Observationally, it has been 
established that the amount of diffuse material is 
nearly ne gligible at th e scal e of individual galax- 
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example, M33, which is ~ 10 11 M^, has been esti- 
mate d to have an /ih l < 1% (|McConnachie et all 
l200l IHood et al.1 1200% . Our own Milky Way is 
thought to have an /mL ~ 1% (e.g., iMorrisor 
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stars from the Sagittarius stream (lLaw et al.l [2 005). 
Our analog M31 has an /ihl ~ 2 — 5% (llrwin et all 
20051; iGuhathakurta et ail 120051: iChapman et all 120061; 
Kalirai et al l 120061) . Recently. IBailin et al.l ([20111 1 



showed that NGC 235, which is of comparable luminos- 
ity to the Milky Way and M31, has /ihl ~ 6%. On group 
scales, ~ 10 13 — 10 14 M , measured /ihl percentages 
tend to be higher than on galaxy scales, though there 
can be considerable variation from group to group (e.g., 



/ihl on the order of a few 
120011: iCastro-Rodriguez et 
/ihl 



iv percen t: iFeldmeier et all 
aTll200l: iFeldmeierl 12006. 



5 - 30%: Da Rocha & Mendes de O liveiral 12005 



lAguerri et ail 120061: iDa Rocha et all 120081 and as high 



as /i hl ~ 45%: IWhite et all 120031: IMcGee fc BaToghl 
2010). For clusters scales, IHL fractions are typically 
much higher. There is still substantia l scatter, but 
ihl v alues range from ~ 10 - 4 0% dMelnick et all 
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Lin fc Mohr l |2004t iMihos et alll2 005; Zib etti et al.l 12005 
Krick etaU \2Q0& iSeigaretall 1200711 As discussed 
in fJ7[ iGonzalez et al.l ([2007D posit that the relevant 
quantity is /ihl+bhg , because it can be difficult to 
disentangle the stellar mass associated with central 
galaxy from the diffuse IHL component. It has been 
found that /ihl+bh g ~ 30% on average on cluster scales 
(IZibetti et aH 120051: IGonzalez et all 120071: iPierini et ail 
2008- lToledo et all 12011 7" 

Our modeling of satellite galaxy stellar mass loss has 
yielded the following principal results and conclusions. 

• Abundance matching with e = and ip = ("tra- 
ditional" abundance matching, akin to using the 
subhalo mass at the time of infall) over-predicts 
the correlation function on small scales (< IMpc) 
for each luminosity threshold sample. 

• For each luminosity threshold, we found both the 
e and ip value that predicted a correlation function 
that matched the data. We found that e and tp were 
decreasing functions of luminosity. This means 
that low-luminosity satellite galaxies are more ef- 
ficiently stripped of stellar material than luminous 
satellites. 

• Our predictions for /ihl from both Model 1 and 
Model 2 are in good agreement with observational 
studies over an enormous range of host halo mass 
scales. For host masses of 10 12 (individual galaxy 
scale), 10 13 — 10 14 (group to small cluster scales), 
and 10 14 — 10 15 (cluster scales), we find mean IHL 
fractions of a few percent, 10 — 20%, and 20 — 25%, 
respectively. These are in accord with previous 
observational studies, though there is substantial 
scatter in the /ihl measurements in the literature. 

• The more physically-motivated Model 2 is con- 
sistent with t he / ihl+bhg measurements of 
IGonzalez et al.l (|2007l ) from small group scales all 
the way through cluster mass systems. The scatter 
in the Model 2 results is comparable to the scatter 
observed in the data. On the contrary, the scatter 
among /ihl+bhg values predicted by Model 1 is 
insufficient to describe the scatter among the ob- 
served systems. This suggests that Model 1, in 
which stellar mass loss occurs in proportion to dark 
matter mass loss, can be rejected by current data. 

We have shown that galaxy clustering can be used as 
a powerful tool to understand how satellite galaxies lose 
stellar mass. We have found the interesting result that 
low-luminosity galaxies lose more stellar mass relative to 
subhalo dark matter mass loss than luminous galaxies. 
Moreover, we were able to predict current IHL obser- 
vations and thus further constrain our stellar mass loss 
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models. These results show that our modeling frame- 
work can mitigate potential problems associated with 
the conventional abundance matching approach. Also, 
our approach is generalizable. It allows for the flexibility 
to make more detailed predictions within this framework 
(for instance, examining the photometric properties of 
the IHL as compared to galaxies). 

In Paper II, we will take advantage of clustering mea- 
surements at high redshifts, allowing us to study the evo- 
lution of satellite galaxy stellar mass loss as a function 
of time. This will enable us to make predictions for the 
build-up of the IHL over cosmic time. A possible ex- 
tension of this program is to also explore stellar mass- 
selected threshold samples directly rather than the lu- 
minosity threshold samples we have used in this paper. 
This is a natural choice for refining this class of stud- 
ies because dynamical models most directly treat stellar 
mass loss. 

A further useful avenue to pursue as a result of this 
work will be to connect detailed theoretical models of 
stellar mass loss in individual galaxies more directly to 



the statistical models of stellar mass loss that can be 
explored with large-scale survey data, such as we have 
done. Models of stellar mass loss that are consistent 
with survey data must also be representative of detailed 
dynamical models of stellar mass loss. It will be very in- 
teresting to develop a set of simple, yet powerful models 
for the build-up and dispersal of the stellar mass in satel- 
lite galaxies that simultaneously describe the evolution 
of galaxy clustering and intrahalo light over cosmic time. 
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